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Abstract 

We present an analysis of the diffuse hard X-ray emission in the core of the 
young massive Galactic cluster Westerlund 1 based on a 48 ks XMM- Newton 
observation. Chandra results for the diffuse X-ray emission have indicated 
a soft thermal component together with a hard component that could be 
either thermal or non-thermal. We seek to resolve this ambiguity regard- 
ing the hard component exploiting the higher sensitivity of XMM-Newton to 
diffuse emission. Our new X-ray spectra from the central (2' radius) diffuse 
emission are found to exhibit He-like Fe 6.7 keV line emission, demonstrat- 
ing that the hard emission in the cluster core is predominantly thermal in 
origin. Potential sources of this hard component are reviewed, namely an 
unresolved Pre-Main Sequence population, a thermalized cluster wind and 
Supernova Remnants interacting with stellar winds. We find that the ther- 
malized cluster wind likely contributes the majority of the hard emission 
with some contribution from the Pre-Main Sequence population. It is un- 
likely that Supernova Remnants are contributing significantly to the Wdl 
diffuse emission at the current epoch. 
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1. Introduction 



Since their launch in 1999, the Chandra and XMM-Newton observa- 
tories have revolutionized the study of X-ray emission from stellar clus- 
ters. Chandra's ACIS and XMM-Newton's EPIC have allowed unprece- 
dented analysis of point sources and diffuse emission in such objects. Of 
particular importance are the observations of extragalactic Super Star Clus- 
ters (SSCs). SSO0 are young (1-10 Myr), massive (10 5 - 10 7 M ) objects 
with extremely dense cores (< 10 5 M pc~ 3 ) and are the predominant sites 
of massive star fo rmation in starburst and interacting gal axies (eg. NGC 



4038/39 and M82, IWhitmore et all . Il999l : iMelo et all 120051 ). However SSCs 



are not limited to these extreme environments, with some found in objects 
such as Blue Compact Dwarfs, non-interacting spirals and Ultra Luminous 
Infrare d Galaxies (eg. Henize 2-10 and PKS 1345+12-C1: iJohnson et al. 



2000l ; lLarsen and Richtlerl . Il999l ; iRodriguez Zaurin et all . 120071 ) . Apart from 
hosting large numbers of massive stars, SSCs also serve to enrich and ener- 
gize the local Interstellar Medium (ISM) through a shocked outflowing cluster 
wind. The cluster wind arises from interacting stellar winds from the massive 
star population and later from SN ejecta. The enrichment of the local ISM 
by the cluster wind can potentially drive further star formation in the region. 
In dwarf starburst galaxies, these winds may be powerful enough to produce 
a galactic outflow, enriching the Intergalactic Medium (IGM) and potentially 
killing further star formation in the galaxy. Hence, SSCs provide not only a 
laboratory for the study of massive stars at various stages of evolution but 
can also provide vital insights into cluster evolution and star formation on 
large scales. Unfortunately, given the distance to many of these SSCs and 
their extremely compact nature it is often impossible to resolve the diffuse 
emission from the point source emission using Chandra or XMM-Newton. 
However, it is possible to resolve the diffuse and point source emission in 
local lower mass analogues. Thus, detailed analysis of such nearer objects 
can provide key insights to the inner workings of SSCs. Westerlund 1 (Wdl) 
is one such cluster, which holds the distinction of being the most massive 
young cluster in the Galaxy. 

Wdl was discovered in the e arly sixties and was initially classified as an 
open cluster (jWesterlundl . I1961I ). The cluster suffers from significant red- 



1 We note that there exists some ambiguity in the literature as to t he classifica t ion o f 
SSCs. We follow here one of the many sets of classification criteria, as in lWhitmore 
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dening [Ay ~ 12.9 mag, iPiatti et all Il998l ) and, because of this, only re- 



cently detailed photometric and sp ectroscopic analyses have been performed 



( Cl ark and Negueruelal. 12001 12004J : iNeeueruela and Clarkl . 120051 : IClark et al. 



20051 ). IClark et all ( 120051 ) f ound a r i ch po pulation of evolved OB stars and, 



using a standard Kroupa (IKroupal . 120011 ) initial mass function (IMF), in- 
ferred a cluster mass of > 10 5 M . This is at the lower limit of the SSC 
mass range and certainly made Wdl the most massive cluster in the Galaxy. 
A more recent deep IR study, h owever, revises this m ass estimate somewhat 
downwards to ~ 4.5 x 10 4 M (IBrandner et all 120081 ) . Although this is still 
bigger than any other known Galactic cluster, it is slightly smaller than ex- 
tragalactic SSCs. The same study also revised previous estimates of age and 
distance to 3. 6±0.7 Myr and 3.55 ±0.17 kpc respectively, which we adopt for 
our analysis. iMuno et all (l2006bf ). henceforth MU06, used Chandra data to 
perform a diffuse emission analysis and found that the emission throughout 
the cluster is dominated by a hard component. However, they were unable to 
identify the nature of this emission due to the absence of hard emission lines 
and discussed both thermal and non-thermal origins for the hard component. 

In this paper we seek to resolve this issue using the XMM-Newton obser- 
vational data, given the telescopes' greater sensitivity to diffuse emission. As 
yet these data have only been used in an analysis of th e well known magnet ar 
CXOU J164710.2-455216 in this cluster flMuno et all l2006ai l2007h . In Sec- 
tion 2 we outline the observational data reduction, before briefly discussing 
the point source analysis in Section 2.1. We follow by presenting the detailed 
diffuse emission analysis in Section 2.2. In Section 3 we discuss our results 
before offering our conclusions in Section 4. 



2. Observations and Analysis 

XMM-Newton observed Wdl on 16 September 2006 for ~48 ks (Obs. 
ID 0404340101, Revolution 1240). The event files were processed using the 
XMM-Newton Science Analysis Software (SAS, Version 7.1.0) meta-tasks 
emproc and epproc. We then filtered the data for good grades in the energy 
band 0.3-10 keV (the energy range at which all 3 of the EPIC instruments 
are most sensitive) and created images for each of the three EPIC cameras, 
namely the PN, MOS1 and MOS2. The PN and MOS2 images were found 
to contain single reflection artifacts which are due to X-rays from a source 
outside the field of view (20'-80' off-axis) reaching the sensitive area of the 
focal plane detectors by single reflection from the rear end of the hyperboloid 
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component of the XMM-Newton mirror shellfl This object was identified 
as the low mass X-ray binary 4U 1642-45 which is located approximately 
20' to the northwest of the observation aimpoint. Images from the three 
EPIC instruments were combined to produce the false colour image shown 
in Figure [TJ 





























1 




Kj..:. - V 




.-IK 







Figure 1: Combined MOS/PN false colour images of Westerlund 1 with red, green and 
blue corresponding to the 0.3-2 keV, 2-4.5 keV and 4.5-10 keV energy bands, respectively. 
North is up, East is left. Left panel: The entire FOV, approximately 30' in diameter. 
Both the Wdl cluster diffuse emission and several point sources are seen at the centre 
with additional sources scattered throughout the FOV. The X-ray binary 4U 1642-45 
reflection artifact is seen in the upper right of the image. Right panel: 5'x5' region 
centered on Westerlund 1 highlighting the cluster diffuse emission. The bright soft source 
to the southwest (seen in red) is the foreground star HD 151018, an 09Iab star. 



2.1. Point Sources 

Point source detection was performed over three standard XMM-Newton 
energy bands (0.5-2 keV, 2-4.5 keV and 4.5-7.5 keV) on the three EPIC 
images using the SAS meta-task edetect- chain. In total, 90 sources with a 



2 See |http : //xmm . esa . int/ external/xmm_user_support/documentat ion/uhb/ 
node23 .html 
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minimum maximum-likelihood detection threshold of 10 were found; 7 of 
these are associated with the reflection and were thus ignored. A further 
8 sources were flagged as spurious due to their positions on or near chip 
gaps and were removed from consideration, leaving 75 source detections in 
the field. By cor relation with the comprehensive Chandra source list in 



Clark et al.l (120081 ) . 4 of our XMM-Newton sources appear to have high mass 



stellar X-ray emitting counterparts in the cluster with a further 8 having 



Muno et al. 




2006b) 


Clark et al. ( 


2008 





2006b) was found to have 



area (within 5' of the cluster centre, 
no counterpart in the source list of 
database and is likely a newly detected flaring PMS star. Table [T] gives our 
detected cluster sources and their corresponding Chandra designations, along 
with spectrally derived source parameters. 

2.2. Diffuse Analysis 

It is obvious from Figure [T]not only that the reflections could contaminate 
the diffuse emission in Wdl but also that they are more prominent in the 
harder energies which are of particular interest to our analysis. To address 
this we considered various analysis techniques including the XMM-Newton 
Extended Source Analysis Software (ESAS) and 'blank sky' background event 
files but found that none could adequately account for the reflection emis- 
sion. Instead we opt for the more traditional method of background extrac- 
tion from regions within the FOV. By defining background regions that are 
as contaminated by the reflection as the cluster, the contribution of the re- 
flection to the cluster spectra can be reduced. We decided against using the 
PN data for the diffuse emission analysis as several detector gaps mask some 
of the Wdl cluster core. Therefore, the following analysis is based on the 
MOS data only. To assess first the diffuse emission in the FOV of the MOS 
cameras we create a non-background subtracted image of the emission using 
ESAS (Version 2]|§. To exclude the point sources from the image we use an 
adapted form of the cheese task to mask not only the po int sources detecte d 
in our XMM-Newton analysis but also those detected by IClark et al. ([2008)0 
but undetected in our XMM-Newton analysis. For the detected sources we 
determined appropriate exclusion regions by assessing the radial brightness 



3 Software available at http://xmm2.esac.esa.int/external/xmm_sw_cal/ 
background/epic_esas.shtml 

4 Elcctronic catalogue accessed via the Vizier service 



5 



Table 1: Westerlund 1 XMM- Newton Point Sources. 



High Mass Stars 



No. 


Chandra Source 
(exo J) 


Opt. ID 


Sp Type 


X-ray Type 
(a) 


MOS Net Counts 
(b) 


Model 

(c) 


N„ 
(10 22 cm -2 ) 


kT 
(keV) 


kT 2 
(keV) 




(io- 


punabs 
X 

12 erg cm -2 s~ 1 ) 

(d) 


T unabs 
(lO 33 erg s _1 ) 
(d) 


1 


164704.1-455039 
164704.1-455031 
164705.1-455041 


W 30 
W 9 
W 27 


O9-B0.5Ia 
sgB[e] 
OB SG 


CWB 
CWB 
RDIS 


2157 


2T 


2.5 


0.7 


2.9 


2.170 




4.37 


6.59 


2 
3 
1 


164705.0-455225 
164708.3-455045 
164710.2-455217 


WR F 
WR A 


WC9d 
WNTb 


CWB 
CWB 
Magnetar 


509 
988 
1708 


IT 
IT 
BB 


0.77i|| 

u ' yu 0.76 


4.72™° 
1.48i;f? 
0-63°;|6 




0.909 
1.967 
1.084 




0.23 
0.89 
0.36 


0.26 
1.34 
0.54 


Pre-Main Sequence Stars 


No. 


Chandra Source 

(exo j) 


Opt. ID 


Sp Type 


X-ray Type 
(a) 


MOS Net Counts 
(b) 


Model 

(c) 


N H 
(10 22 cm- 2 ) 


kT 
(keV) 


kT 2 
(keV) 


x7" 


(io- 


runabs 
X 

12 erg cm -2 s~ 1 ) 

(d) 


t unabs 

^x 

(lO 33 erg s- 1 ) 
(d) 



0.959 0.04 0.05 

1.992 0.01 0.02 

1.057 0.06 0.09 

0.954 0.12 0.19 



In this table the Chandra designation, optical ID and spectral types are adopted from lClark et af . (2008) and references therein. Notice that Source 
1 encompasses W 30, W 9 and W 27 which are unresolved by XMM-Newton. No counterpart for Source 7 was found in the analysis of lCIark et al.l 
d200STl or in the SIMBAD database. 

(a) - X-ray type adopted from IClark et all 1120081 ). PMS = Pre-Main Sequence, CWB = Colliding Wind Binary, RDIS = Radiatively Driven 
Instability Shocks in stellar winds, (b) - Combined MOS net counts. PN counts omitted due to some sources lying on detector gaps, (c) - Best fit 
X-ray spectral model each of which was convolved with an absorbing hydrogen column. Only those sources with > 100 combined MOS and PN net 
counts were spectrally analysed. All spectral fits were performed in XSPEC 12.3.1. IT = one temperature thermal plasma, 2T = two temperature 
thermal plasma, BB = black body. The spectrum of Source 2 is approximated by a 2T thermal plasma model however the true spectrum is likely 
described by a more complex model given it is an amalgam of several unresolved bright sources. In addition, Sources 3 and 4 are likely described 
by a 2T thermal plasma model however the soft components are poorly constrained and thus only the IT approximation is given, (d) - Absorption 
corrected X-ray fluxes and luminosities quoted in the 0.5-8 keV energy range. 



5 
6 
7 
8 
9 
10 
11 
12 
13 



164640.8-454834 
164648.8-455307 

164652.6-455357 
164703.2-455157 
164712.8-455435 

164713.6- 454857 

164718.7- 454758 
164720.1-455138 



PMS 
PMS 
PMS 
PMS 
PMS 
PMS 
PMS 
PMS 
PMS 



Flare 
Flare 
Flare 
Flare 
Flare 
Flare 
Flare 
Flare 
Flare 



34 
39 
62 
210 
397 
10 
62 
66 
19 



IT 
IT 
IT 
IT 



0-76J;! 
0.5l0;| 



°- 86 25 

4.07^1 
2.491:1° 



profile of each source in two energy bands (namely 2-4.5 keV and 4.5-10 keV, 
since we will focus on the harder energies in the analysis below). In each 
energy band the extent of the exclusion regions is set by the source bright- 
ness profiles approaching local background levels and the largest of these two 
areas is adopted as the overall source exclusion region. For each of the XMM- 
Newton sources in the Wdl cluster these regions are ~30" in diameter. This 
is not a perfect solution however since incidental photons from extreme edges 
of the wings of the source PSF may not be masked. We estimate from the 
background subtracted source brightness profiles that no more that 10% of 
the photon flux could have been missed by these exclusions regions. No such 
method can be employed to mask the undetected sources due to the simple 
unavailability of the source brightness profiles. Because of this, 10" diameter 
regions were used to mask as much as possible any detected photons from 
these sources (for which the greatest probability for these photons is to be 
found in the inner PSF core) without sacrificing large numbers of diffuse 
source photons. The resulting image is shown in Figure [2j 

Figure [2] clearly shows that the diffuse emission in the FOV is dominated 
by the reflection. Extracting background spectra from within the FOV would 
mean that any contributing sources of background photons will be contained 
in these spectra, including the reflection itself. Hence, we first assess each 
contaminant as some of these can vary across the FOV. The ESAS man- 
ual identifies 4 sources of background contamination. These are Solar Wind 
Charge Exchange (SWCE) contamination, instrumental fluorescence lines, 
residual soft proton contamination and the cosmic background. In the case 
of this observation, the SWCE and one instrumental fluorescence line can 
largely be ignored as they only significantly contaminate spectra below 1.5 
keV (given the absorbing hydrogen column in the direction of Wdl, any 
emission below this energy is likely foreground emission). The second fluo- 
rescence line (E = 1.75 keV) may contaminate the spectra especially in larger 
extraction regions. Soft proton contamination blights almost all observations. 
Periods of high contamination are screened out in the data reduction process 
but residual contamination can still affect the observational data. However, 
it is assumed that the level of contamination is constant across the FOV and 
hence should be contained in the backgrounds. Similarly it is assumed that 
the cosmic background is constant across the FOV and is also contained in 
the backgrounds. The final source of contamination, unique to this obser- 
vation, is the 4U 1642-45 reflection. This can be compensated for by the 
selection of appropriate background regions that are equally contaminated 
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Figure 2: Smoothed, non-background subtracted, combined MOS1/MOS2 image in which 
sources are masked with the CCDs most affected by the reflection omitted. Again the 
FOV is approximately 30' in diameter. The green circle indicates the region chosen for 
background extraction. This source free region lies on the same CCD as the cluster core in 
both the MOS1 and MOS2 images, is outside the total cluster region and is approximately 
as contaminated by the reflection as the cluster core. A zoomed image of the centre region 
in this figure can be found in Figure [7J 



as the cluster, determined using Figure HJ 

We adopt the cluster core coordinates determined by MU06 as a = 16 
47 04.3, 5q = -45 50 59 and extract spectra from the central 1' radius region 
and three annuli extending out from the core (l'-2', 2'-3.5', 3.5'-5') with the 
XMM-Newton and Chandra point sources masked. As any X-ray photons 
detected below 1.5 keV are likely foreground emission and because of the two 
instrumental fluorescence lines at 1.49 keV and 1.75 keV, we have restricted 
our spectral analysis to the 2-8 keV energy band which conforms to the 
standard Chandra hard band as used by MU06. The MOS 1 and MOS 2 
spectra and ancillary files for each annulus were combined^] and the resulting 
spectra were adaptively binned so that each bin has a S/N of 3. MU06 found 
that the cluster diffuse emission spectra are well fit with either an absorbed 



5 Following the procedure found at |http: //xmm.esa. int/sas/current/ 
|documentation/threads/epic_merging.shtml| 
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two-temperature thermal plasma model (the harder thermal component with 
sub-solar abundance to explain the lack of hard emission lines) or an absorbed 
thermal plasma plus power law model. Hence, the combined MOS spectra 
were fit in XSPEC0 with these models. The abundance parameter of the cool 
thermal component in the fits was poorly constrained due to the majority 
of this emission falling below 1.5 keV, however it was fixed at 2 Z to be 
consistent with MU06. The fit results are in general agreement with those of 
MU06. The outer annulus exhibits only a hard component adequately fit by 
either an absorbed thermal plasma or absorbed power law model. The best- 
fit absorbed two temperature thermal plasma and absorbed thermal plasma 
plus power law models yield very similar best-fit statistics in the inner annuli 
due to the absence of hard emission lines, as also found by MU06. This is 
somewhat surprising as, for a cluster such as Wdl with a large, centrally 
located WR population, the hard diffuse emission is e xpected to be th ermal 
in origin due mostly to a thermalized cluster wind ( lOskinoval . 120051 ). To 
assess this further we extracted and combined MOS spectra from the inner 
T radius region which again was adaptively binned so that each bin has a 
S/N of 3. The resulting combined MOS spectrum is shown in Figure El 



With this choice of extraction region, the He-like Fe 6.7 keV line, which is 
used as a diagnostic for thermal diffuse hard emission, is now seen in the com- 
bined MOS spectrum. A straightforward summation of the line flux above 
continuum, combined with the counts statistics, yields a line significance of 
6.2 a. We fit the 2-8 keV combined core MOS spectrum in XSPEC with 
various combinations of thermal and non-thermal models and find that the 
data are best fit with an absorbed two temperature thermal plasma model 
as this can account for the He-like Fe 6.7 keV line and the softer emission 
lines below 3 keV due to He-like S and He-like Si. The results of the core two 
temperature thermal plasma and, for comparison, the thermal plasma plus 
power law fits are given in Table [21 together with values for the diffuse X-ray 
surface brightness corresponding to these fits. The fact that the two temper- 
ature thermal plasma model is a better fit as only this can account for the 6.7 
keV emission line demonstrates that, at least in the inner 2' radius region, 
the hard component is predominantly thermal in origin. We applied a similar 



6 XSPEC Version 12.3.1 was used for all spectral fits in this analysis. 
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Table 2: Diffuse Spectral Fits 



Region 


(lO 22 cm" 2 ) 


(kcv) 


Two temperature thermal plas 
(Z/Z )jI fcT 2 (Z/Z Q ) 2 

(kcv) 


;ma 


I 10 erg cm 


6 


r unabs 6 
(lO 33 crg^ 1 ) 


(lO 3 


T unabs 
^X,SB 

3 erg s -1 


C 

PC" 2 ) 


<2' 


^■ uo 1.88 




2 


o n7 3 ' 67 

°- u '2.69 


U - DZ 0.40 


0.971 


1.71 




2.56 




0.27 




Region 


(lO 22 cm- 2 ) 


(kcv) 


Thermal plasma plus power law 

(z/z )i Q r - x 2 /v 


(l0 — 12 erg cm" 


6 

-».-o 


r unabs b 
^X 

(lO 33 ergs- 1 ) 


(io 3 


t unabs 
-^X,SB 

3 erg s- 1 


c 

PC" 2 ) 


<2' 


9 fl7 2 - 32 
/ - u ' 1.81 


0.81°;72 


2 


2-43l;l? 




1.142 


1.50 




2.27 




0.24 





a (Z/Z Q )i for both models is fixed at 2. 

6 Absorption corrected X-ray fluxes and luminosities quoted in the 2-8 keV energy range, consistent with the analysis of 
MU06. 



c Diffuse emission surface brightness in the 2-8 keV energy range 



2 



5 

Energy (keV) 



Figure 3: Combined MOS inner 2' diffuse emission spectrum binned with S/N=3 and 
fit with an absorbed two temperature thermal plasma model. The blue and red lines 
indicate the individual thermal components of the model, representing the hard and soft 
components respectively. 



treatment to both the outer extraction annuli, however this failed to reveal 
any line emission in the hard continuum so it is still debatable as to whether 
thermal or non-thermal processes are responsible for the hard emission in 
these regions. Diffuse non-thermal X-ray emission has been found for only 



a few star formation regions in the Galaxy and in the LMC (IMaddox et al 



20091 . and references therein). As a successful fit was obtained for the Wdl 



diffuse emission including the Fe 6.7 keV (Table |2]), this leaves little room 
for any additional non-thermal contribution. Nevertheless, we have run sev- 
eral further, explorative spectral fitting attempts with XSPEC for models 
with thermal plus non-thermal components. However, even when fixing non- 
thermal fitting parameters to reasonable choices, no convergence of the fits 
could be obtained. Thus, no evidence for an additional diffuse non-thermal 
emission component was found for Wdl. 
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3. Discussion 

3.1. Origin of 6.7 keV Emission Line 

To confirm that the emission line at 6.7 keV is a feature of the cluster 
diffuse emission we must rule out other potential sources for this line, such 
as point source contamination, the reflection, cosmological background and 
those other contributions mentioned above. 

3.1.1. Point Source Contamination 

Although exclusion regions were defined to mask the point sources it is 
possible that these masks did not completely exclude point source photons 
at the extreme edges of the wings of their PSFs. To verify that the emission 
line is not the result of point source contamination, combined point source 
spectra for those sources that lie within the 2' radius circle centered on the 
cluster core were extracted for each MOS camera. The resulting spectra were 
combined and this spectrum adaptively binned so that each bin has a S/N 
of 3, shown in Figure HI 




2 5 
Energy (keV) 



Figure 4: Combined MOS spectrum of point sources within a 2' region centered on the 
cluster core. 

Inspection of the 6.7 keV region of the combined MOS point source spec- 
trum in Figure S] shows line emission centered at ~ 6.7 keV with a significance 
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of 7.3 a. To assess the strength of this line in comparison to the 6.7 keV line 
in the diffuse emission spectrum (which excludes the point sources) we de- 
termine the X-ray flux in the 6.6 keV - 6.8 keV range for the combined MOS 
diffuse spectrum and the combined MOS point source spectrum. We find 
that the combined MOS point source spectral line flux of 5.73 x 10~ 14 erg 
cm -2 s" 1 is just over half the strength of the combined MOS diffuse emission 
line flux of 1.02 x 10 -13 erg cm -2 s _1 . If the emission line in the diffuse 
spectrum were a result of point source contamination we would expect the 
line flux in the point source spectrum to be significantly larger than that in 
the diffuse spectrum given that the majority of these photons was in fact 
masked out when creating the diffuse spectrum. The exclusion regions for 
the sources detected in the XMM-Newton observation were defined based on 
their brightness profiles in the 2-4.5 keV and 4.5-10 keV energy bands (see 
Section 12. 2j) with the largest determined region from these bands set as the 
overall exclusion region for the respective sources. The exclusion regions were 
invariably determined by the 2-4.5 keV brightness profiles given the higher 
number of photons in this energy band in comparison to the 4.5-10 keV band 
for each of the sources. In addition it was determined that no more than 
10% of the source photon flux is missed by these extraction regions. Since in 
this discussion we are focussing on the 6.6-6.8 keV energy range and given 
that for sources with off-axis angle < 2' in EPIC MOS observations (which 
is our Wdl core radius) the fractional encircled energy radii decrease with 
increasing photon energ}0, the exclusion radii for the sources should be more 
than sufficient to mask >90% of the source photons in the 6.6-6.8 keV energy 
range. As such we expect the point source contamination due to photons from 
the extreme wings of their PSFs to contribute < 5% of the observed diffuse 
emission line flux. Additionally, the point source extraction regions occupy 
~25% of the cluster core area which is one third of the diffuse emission ex- 
traction region. As such one third of the combined diffuse emission flux in 
the 6.6 keV - 6.8 keV range is likely contained in the combined point source 
spectral line flux, meaning the flux due to point sources only may be as low 
as 2.33 x 10 -14 erg cm" 2 s" 1 and contribute <<5% of the observed diffuse 
line flux. For these reasons we can rule out point source contamination. 



See |http : //xmm . esa . int/ external/xmm_user_support/do cumentat ion/uhb/ 
|node!8 .htmll 
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3.1.2. Reflection 

To reduce the contribution of the reflection to the diffuse spectra we se- 
lected background regions as much contaminated by the reflection as the 
cluster. However, this is not a perfect solution and contamination by the re- 
flection has to be investigated. To assess if contamination by the reflection is 
the source of the emission line, we extract and combine spectra for those re- 
gions most affected by the reflection (but without the reflection hyperbolas). 
We create a new combined MOS background spectrum from regions devoid 
of point sources on the opposite side of the FOV where the contamination 
from the reflection is at a minimum (i.e. - mainly the regular background 
contaminants present) and subtract this from the combined MOS reflection 
spectrum. The resulting spectrum was again binned so that each bin has a 
S/N of 3, shown in Figure [5j 




2 5 
Energy (keV) 



Figure 5: The combined MOS background-subtracted reflection spectrum, extracted from 
those regions most contaminated by the strong NW out-of-FOV source, is shown in black. 
The background spectrum, shown in red, is extracted from regions devoid of point sources 
on the opposite (SE) side of the FOV, where the contamination from the reflection is at a 
minimum. The background spectrum was binned so that each bin contains 30 counts. 

Figure [5] shows no significant emission at or around 6.7 keV so we can 
safely assume that contamination by the reflection is not responsible for the 
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line. 



3.1.3. Instrumental Background 

The instrumental background comprises the soft proton background, in- 
strumental fluorescence lines and particle background. From the ESAS method, 
it is seen that the soft proton contamination is well modeled as a power law 
(or broken power law). Hence soft proton contamination cannot be respon- 
sible for the emission line. Similarly the instrumental fluorescence lines are 
modeled as Gaussian peaks at 1.49 keV and 1.75 keV and thus are not the 
source of the line at 6.7 keV. Moreover, the p article background contribu tion 



is simply too low to cause the emission line ( jKuntz and Snowdenl . 120081 ) . 



3.1. 4- Cosmic Background 

When extracting the backgrounds we assume that the cosmic background 
emission is constant across the FOV and thus contained in the backgrounds. 
To assess the contribution of the cosmic background we use the combined 
MOS background extracted for use with the reflection spectrum in Section 
I3.1.2[ indicated by the red spectrum in Figure [5j We see no significant 
emission lines at or around the 6.7 keV region so we can assume that the 
diffuse emission line is not from the cosmic background. 

3.2. Sources of Hard Emission 

Having demonstrated that the 6.7 keV line is a feature of the cluster 
diffuse emission and hence that the diffuse emission within 2' of the cluster 
core is mostly thermal in origin, we now address possible sources for this 
component. There are three potential sources of thermal emission in the 
core of a cluster as large as and with the age of Wdl. These are unresolved 
PMS stars, a thermalized cluster wind and Supernova Remnants (SNRs). 

3.2.1. Unresolved Pre-Main Sequence Stars 



In their Chandra point source analysis, Oark et al.l ( 120081 ) identified not 
only the X-ray emitting high mass population but also a number of PMS 
stars down to a limiting luminosity of ~ 10 31 erg s -1 . These sources repre- 
sent the most luminous members of the PMS population and were masked in 
the spectral analysis above. However, the remaining unresolved PMS pop- 
ulation could contribute significantly to the diffuse emission in the cluster. 
To obtain an estimate for the contribution of unresolved PMS population, a 



15 



similar method to t hat adopted bvlGetman et al.l ( 120061 ) ; IWang et al.l ( 120071 ) ; 

Broos et all (l2007t ) : lEzoe et all ( l2006l )7 and others was used. 

We constructed a PMS source 0.5-8 keV photo n flux distribution (Fx, phot; in 



units of phot cm 2 s 1 ) from the catalogue of Park et al.l ( 120081 ) and com- 



pared it to that of the Orion Nebular Cl uster (ONC) using d ata from the 



Chandra Orion Ultradeep Project (COUP, iGetman et al.l . 120051 ) . We assume 



that the Wdl cluster and ONC have the same Initial Mass Function (IMF) 
and thus X-ray flux distribution, differing only in the size of their underly- 
ing populations. However, prior to this, several corrections were applied to 
COUP photon fluxes to account for the increased distance, foreground ab- 
sorption and age of Wdl. To account for the difference in distance the COUP 
fluxes were simply scaled according to the ratio of the cluster distances. Since 
the ONC and Wdl are s ubject to quite differen t absorbing hydrogen column 



densities (3 x 10 21 cm 2 (IFeigelson et al.l . 120051 ) and ~ 2 x 10 22 cm 2 , respec- 



tively) we must correct the COUP photon fluxes for the increased absorption 
to Wdl. We do this by assuming that the composite spectrum of the Wdl 
PMS population is the same as that of the COUP, namely a t wo-temperature 
therm al plasma with kTi _= 0.5 keV and kT 2 = 3.3 keV (IFeigelson et al 



2005 ). Using the PIMMS0 tool we determine an absorption correction for 
the composite spectrum and apply it to the COUP photon flux distribution. 
Finally, since the X-ray luminos ity of PMS stars has been foun d to decrease 
with time r as L x oc r~ - 3 ( jPreibisch and Feigelsonl . |2005|), the COUP 
fluxes were adjusted to account for the decreased X-ray luminosity of the 
older PMS population in Wdl. At this stage, the COUP photon flux dis- 
tribution has been adjusted so that the ONC population effectively has the 
same age, distance and absorption as Wdl. The photon flux distributions 
can now be reliably compared, shown in Figure El 

Figure M indicates that the higher end of the flux distribution of each of the 
clusters is very similar, differing only in the number of sources per bin, which 
shows that the earlier assumption of a similar underlying IMF, and hence X- 
ray flux distribution, was reasonable. The break in the Wdl distribution 
just below log F XiP hot ~ — 6 phot cm -2 s -1 is due to the incompleteness of 
the stellar population because of observational constraints. To determine the 
contribution of the undetected population we scale the COUP distribution 
to that of Wdl to determine the overall photon flux of the entire population 



See Ihttp : / /heasarc .nasa.gov/docs/software/tools/pimms .html 
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Figure 6: Photon flux distributions for the Wdl and ONC. The ONC distribution has 
been corrected for the increased distance, foreground absorption and age of Wdl. The 
errorbars indicate the la Poisson error. 



and subtract that of the detected population, leaving a photon flux due to 
the unresolved sources of 4 x 10 4 phot cm" 2 s -1 . We can convert this to a 
luminosity if we consider that the overall absorption corrected X-ray lumi- 
nosity of the COUP population responsible for the photon flux distribution 
is known (overall 2-8 keV L^ ahs ~ 1.2 x 10 33 erg s _1 , corrected for the age 
of Wdl). Scaling this to the unresolved Wdl population yields a 2-8 keV 
Ly nabs of ~ 1 3 x 10 33 erg s _1 . However, not all of the PMS population of 



Clark et al.l ( 120081 ) is located in the core of the cluster but rather spread 



through the entire cluster area. Since the vast majority is within 3 pc of the 
cluster centre, we assume that all of the determined unresolved PMS X-ray 
emission comes from this region, resulting in a surface brightness value of 
~ 4 x 10 31 erg s~ x pc~ 2 . This is approximately 15% of the surface brightness 
of the diffuse emission in the core (see Table |2]), meaning an unresolved pop- 
ulation can only account for a small fraction of the diffuse emission in the 
core. 
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3.2.2. Cluster Wind 



A cluster as large as Well is expected, and is indeed found by I Clark et al 



( 120051 ). to contain many massive stars in the cluster core. These massive stars 
are the source of large amounts of energy and mass being ejected into the 
cluster volume via stellar winds. The winds collide and thermalize, filling the 
cluster core volume with a hot, shocked plasma. After some time, the outflow 
from these th ermalized winds be comes stationary and a steady state cluster 



wind ensues ( jCanto et al.l . 120001 and references therein). The temperature 
of this hot, diffuse plasma throughout the core is sufficiently high to radiate 
at X-ray energies and hence is potentially responsib le for the diffuse h ard 
emission in Wdl. MU06 make use of the equations of ICanto et al.l (120001 ) for 
the central hydrogen density and temperature, namely: 



n 



CHI" 



0.1N 



M 



1O" 5 M yr" 



10 3 km s 



-i 



Rc 
pc 



^ (1) 



To 
K 



1.55 x 10' 



10 3 km s 



(2) 



where N is the number of stars contributing to the thermalized cluster wind, 
M is the average mass loss rate per star, v w is a weighted average wind 
velocity of the stars and Rc is the radius of the region containing the stars. To 
enable direct comparison of results we apply a similar treatment as performed 
by MU06. In the analysis of MU06 only the WR stars are considered as the 
contributors to the cluster wind and from Equation 1 an emission measure 
(Kem = ^ 7! 'Rc n o) * s calculated before using a standard thermal plasma model 
in XSPEC to extract fluxes. Our treatment differs slightly from that of 
MU06 in two ways. First, we set the value of Rc at 2 pc (~2' at 3.55 
kpc) which is our core extraction reg ion in the analysis above. This is an 
acceptable value as [Canto et al.l (l2000f ) state that the value of Rc needs only 
be approximate to the distance from the cluster centre to the outermost star, 
i.e. - all the stars considered to contribute to the cluster wind are inside Rc 
(see also Figure Cj). Second, rather than assign typical values for M and 
v w for WRs i n general, we use the known spectral types of the WRs in the 
cluster core dCrowther et al.l. 120061) along with general physical and wind 
properties of WRs (jCrowtherl . 120071 ) to estimate the more accurate mean 



values of M = 1.4 x 10" 5 M Q yr" 1 and v w = 1320 km s" 1 for the 21 WRs 
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within 2 pc of the cluster centre. Inputting these values into Equations 1 
and 2 yields n = 0.5 cm -3 and kT = 3.7 keV. We then c alculate K F m 



and f eed the appropriate values into the APEC spectral model ([Smith et al 



20011 ) in XSPEC which gives an unabsorbed X-ray luminosity (Lx nabs ) in the 
2-8 keV energy range of ~ 2 x 10 33 erg s _1 . Extracting the observed high 
temperature value from our two temperature fit yields an L x nabs of 1.7 x 10 33 
erg s _1 . Hence, we find the predicted value is in excellent agreement with 
the observed value. Thus, a thermalized cluster wind can account for the 
hard thermal emission in the core of Wdl. Obviously, this remains true also 
if the estimated contribution by unresolved PMS objects (Section 13.2 . 1 [) is 
subtracted from the observed hard emission. 

This result differs from that of MU06, who found that the hard emission 
in the cluster core is approximately half of their predicted value. The dis- 
crepancy results mainly from the setting of Rq = 4 pc (~3 ! at their adopted 
distance of 5 kpc) in their calculations. Our smaller Rq value reduces the 
emission measure and hence the derived L x nabs . We note further that the 
morphology of the Wdl diffuse X-ray emission (see Figure [7]) resembles the 
distribution of the massive stars in the cluster, in particular in showing a 
noticeable extension to the SE. When the less massive stars are considered, 
such as those down to O.8M (Brandner et al. 2008; see also their Figure 
1), the cluster shows a smoother, roughly N-S elongated distribution. Also 
these qualitative considerations link the diffuse emission with the massive 
sta rs (that exhibit winds ) rathe r than the low-ma ss (PMS) ob j ects. 



Stevens and Hartwelll ( 120031 ) used the model of ICanto et al.l ( 120001 ) to pre- 



dict the properties of the cluster wind for some Galactic and Magellanic Cloud 
clusters and compared them to observation. They found that, in all cases 
bar one, the predicted kT values are much larger than observed and that 
the predicted X-ray luminosities are much less than observed. Our results 
are somewhat at odds with those of Stevens h Hartwell in that our predicted 
values of kTo and L x nabs are quite close to those observed. The difference in 
kT can be explained by our adopted v w value. For the determination of kT 
we are only considering the WRs in the cluster (given that the early O stars 
have already evolved off the main sequence), which have lower wind velocities 
than the early O stars in the clusters studied by Stevens & Hartwell, thus 
keeping kTo down. A possible explanation for their difference in L x nabs , and 
also kTo values, is that in all bar the same one cluster, Stevens & Hartwell 
fit the diffuse spectra with a single temperature model. This may underes- 
timate any cool thermal component that is present in the diffuse spectra of 
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Figure 7: Distribution of massive stars in the cluster core in relation to the diffuse emission. 
Wolf-Rayets are indicated by the filled green circles whereas main sequence and evolved 
OB stars are indicated by filled blue circles. The large circle indicates the 2" core radius 
and the black areas indicate areas where point sources were masked. It is clear that the 
diffuse emission resembles the distribution of the massive stars in the cluster core. 



the clusters and the observed kT value will be an amalgam of the hot and 
cool component kT values (in fact, results for kT and Lx nabs from single tem- 
perature fits will be coupled). If a cool component has been underestimated, 
this will dramatically affect the observationally derived L^ ahs given that the 
unabsorbed cool component will greatly add to the overall Lx nabs . Since we 
allow for a cool thermal component in Wdl and restrict our analysis to the 
harder thermal component only, our results are close to those predicted. One 
must be aware however that the model used in our analysis does not incor- 
porate the thermalization efficiency and mass loading of the cluster wind 
described by Stevens & Hartwell. In practice, our adopted model assumes no 
mass loading and a thermalization efficiency of 1 (i.e. - no radiative losses in 
the conversion of the stellar wind energies to the cluster wind). This may be 
simplistic and we note here that a change in either parameter would serve to 
increase the predicted overall X-ray luminosity (i.e. - not restricted to the 
hard energy band) and reduce the kT value. 
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3.2.3. Supernova Remnants 

SNRs emit X-rays through thermal and/or non-thermal processes. Since 
we have identified the hard emission in the core of Wdl as predominantly 
thermal, we address only the thermal emission mechanisms (i.e. - the SNR 
interaction with its surroundings). The occurrence of SuperNov ae (SNae) in 



Wdl has been considered in previous analyses of Wdl (MU06; Park et al. 



20081 ; iBrandner et all 12008k etc.), which use the stellar population to ex- 
trapolate the cluster IMF to higher masses and determine that the cluster 
initially contained ~ 10 2 stars with M > 50M Q . Given the age of Wdl we can 
assume that all these stars have already been lost to SNae. However, apart 
from the magnetar, no evidence of post-SN objects (SNRs, com pact objects 



or X-r ay Binaries) was found in the analyses of either MU06 or [Clark et al. 
( 120081 ). Possible reasons for the absence of the discrete objects are beyond 
the scope of this analysis, instead we only address the potential contribution 
of any SNR to the diffuse emission in the cluster core. SNRs are expected 
to emit X-rays when the shock front interacts with the surrounding ISM, 
but in the Wdl region winds from the massive stars have cleared away the 
ISM so we cannot expect to observe emission from this process (MU06 and 
references therein). 

However, SNRs can also interact with the stellar winds of nearby massive 
stars in the cluster. An illustrative exam ple is the LBV/W R binary system 
HD 5980 in the Small Magellanic Cloud (INaze et all . 120021 ). the stellar wind 
of which is likely interacting with the ejecta from the nearby object SNR 
0057-7226. INaze et al.l (120021 ) performed an analysis of the diffuse emission 
around the HD 5980 system finding several X-ray bright filaments in the dif- 
fuse structure and obtaining a 0.3-10 keV L™ abs of ~ 10 35 erg s -1 , result s that 
both were subsequent l y deri ved theoretically by I Velazquez et al.l ( 120031 ). Ad- 
ditionally INaze et al.l (120021 ) found that the diffuse emission associated with 
the object is best fit by an absorbed thermal plasma with kT = 0.66 keV. 
This value is much lower than the hard component plasma temperature of 
the diffuse emission in Wdl, already suggesting that SNRs are not primarily 
responsible for the hard diffuse emission in the core since this plasma tem- 
perature cannot account for the obser ved 6.7 keV emiss ion line. However, in 
view of the large luminosity found bv lNaze et ail ( 120021 ) for the diffuse emis- 
sion around HD 5980, it is possible that the high energy tail of a SNR/stellar 
wind interaction spectrum may contribute to the observed hard diffuse emis- 
sion in Wdl. To assess this we consider the cool thermal component of the 
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Wdl core diffuse spectrum which has a derived plasma temperature of 0.68 
keV, similar to that of the HD 5980 diffuse emission. From the results of the 
spectral fits in Table [2] we estimate that the soft component contributes about 
one third of the 2-8 keV Lx nabs in the core. Or, if the estimated hard emission 
contribution by unresolved PMS objects of 0.46 — 1 x 10 33 erg s _1 (Section 
I3.2.ip is subtracted from the observed hard diffuse emission, about 40-60% 
of the remaining hard emission could be contributed by the soft component. 
Thus even if SNRs were responsible for all of the observed soft emission, they 
would only account for between one third to one half of the observed hard 
emission in the core but crucially could not account for the 6.7 keV emission 
line. 

In itself, any SNR contribution to diffuse emission is of cours e dependent 



on the recent occurrence of a SN event in or near the core. Muno et al 



( I2006bl ) determine the SN rate in Wdl to be once every 7,000-13,000 yr. 
If we assume that a SN event occurs at the very centre of the cluster core 
and that the SN ejecta travel with a velocities of a few 10 3 km s _1 , then 
the ejecta escape the core region after several 10 2 yr (and the cluster itself 
after a few 10 3 yr). This leaves a significant amount of time (typically over 
90% of the time between successive SNae) in which the core is free from the 
effects of SNRs. Also this time scale argument suggests that it is unlikely 
that SNRs are contributing significantly to the Wdl diffuse emission at the 
current epoch. 

4. Conclusions 

The above analysis and discussion of the until now unexploited XMM- 
Newton diffuse emission data for Wdl have demonstrated a Fe 6.7 keV emis- 
sion line, indicating that the hard X-ray component in the inner 2' radius 
region of the cluster is predominantly thermal in origin. The most likely 
explanation for this diffuse component is a thermalized cluster wind. An es- 
timated value for the 2-8 keV X-ray luminosity produced by a cluster wind in 
Wdl (2 x 10 33 erg s _1 ) is close to the observationally determined value for this 
luminosity (1.7 x 10 33 erg s _1 ). The conclusion that the cluster wind is the 
li kely cause of the diffuse emission is also in line with the model predictions 



of lOskinoval (120051 ) . 

While the unresolved PMS objects are less likely the main cause of the 
Wdl diffuse emission, they could nevertheless still make a contribution to the 
diffuse component. We determined the unresolved PMS source contribution 
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by comparing the observed PMS source flux distribution of Wdl to that of the 
ONC. We found that the derived surface brightness value for the unresolved 
source population of ~ 4 x 10 31 erg s" 1 pc -2 can only account for ~ 15% of 
the observed value of ~ 3 x 10 32 erg s _1 pc -2 . 

SNRs interacting with stellar winds are likely to be too soft X-ray emitters 
than is required to produce the 6.7 keV emission line. Individual SNRs have 
not been identified in Wdl and, from time scale arguments, are a priori not 
likely to be present. 
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